Renormalization of 1 Sq NN scattering amplitude in Effective Field Theory 
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Cutoff regularized subleading order So NN potential of effective field theory (EFT) is 
iterated using Lippmann-Schwinger equation. It is shown that the scattering amplitudes 
calculated in cutoff and subtractively renormalized EFT are equal up to the accuracy 
of performed calculations. Non-perturbative renormalization, where part of divergences 
are absorbed into two contact interaction coupling constants with subsequent removal 
6JT), of regularization is also performed. Cutoff and dimensional regularizations both lead to 

3 ■ finite but different results within this scheme. 
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(N : 

There has been significant development in effective field theory (EFT) approach to problems of multi 
nucleon systems originated by Weinberg's papers PJ, 0- (for recent developments see refs. [|J, [[|). In 

■ Weinberg's approach power counting is formulated for effective potential and the scattering amplitude 
is calculated by substituting this effective potential into Lippman-Schwinger (or Schrodinger) equation. 
In rcf. H it has been argued that Weinberg's approach suffers from formal inconsistences. While this 
conclusion has been questioned @] , @ a new explicitly consistent power counting (KSW power counting) 
has been suggested in ref. Q . In this approach pions are included perturbatively and it allows to perform 

■ all calculations analytically which can be considered as a considerable advantage. Unfortunately it has 
been realized that perturbation theory series do not converge in this approach Q , ^ , . An attempt 

• to combine the advantages of Weinberg's and KSW schemes has been made recently in ref. |pT| . In that 
paper it has been confirmed that Weinberg's power counting is formally inconsistent in the channel 
NN scattering problem. In particular, divergences which appear in leading order, can not be absorbed 
by leading order counter-terms. Different approach to the problem of non-perturbative renormalization 
of Lippmann-Schwinger equation has been suggested recently in ref. p2[ . 

In present letter, which is the continuation of the work started in refs. fl3j| , [fbif , we follow Weinberg's 
approach and iterate sub-leading order 1 So NN potential in the framework of Lippman-Schwinger equa- 
tion. We note that the above mentioned inconsistency of Weinberg's power counting manifests itself only 
when it is applied to unrenormalized diagrams and counter-terms separately. We advocate the point 
of view that the power counting should be applied to renormalized diagrams @, Q], [jTsf , ijlGf . Let us 
emphasize that we do not attempt to resolve the above mentioned formal inconsistency of Weinberg's 
power counting, we believe that it can not be resolved within this scheme. We share the point of view 
that this problem is irrelevant: if one applies the power counting to renormalized diagrams then there 
are no formal inconsistencies, renormalized diagrams are of the order which has been assigned by power 
counting and the approach is successful in studies of processes in few nucleon systems [|l7| , [p| , [|l9| , |Q . 
This scheme can be implemented using BPHZ renormalization procedure. In BPHZ approach no counter- 
terms are included into the Lagrangian but instead the Feynman rules are accompanied by prescriptions 
for subtracting loop diagrams pl| . The subtracted diagrams depend on renormalization point (s). This 
dependence is cancelled by the renormalization point dependence of renormaliscd coupling constants so 
that the physical quantities are independent of the choice of subtraction point (s). The above mentioned 
problem of inconsistency that higher order counter terms contribute into the renormalization of low order 
diagrams translates now into the following problem: low order diagrams give large contributions into beta 
functions of coupling constants of higher order (counter) terms. There have been concerns that even if 
for some value of the renormalization point (s) the higher order coupling constants are tuned to be small, 
the slight change of this renormalization point (s) would result in large values for coupling constants and 
therefore power counting would be badly violated. Note that these concerns arc based on pcrturbative 
estimations of beta-functions in the area where perturbation theory is not actually applicable. Hence 
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there is no reason to expact that the assumption of naturalness which is an ingredient of Weinberg's 
power counting is violated. To make any reliable conclusions about renormalization group behaviour of 
renormalized coupling constants non-perturbative analysis of this behaviour are required. We are not 
aware of any such analysis of renormalized coupling constants in EFT with explicitly included pions. On 
the other hand the considerable success of cutoff EFT, which is based on Weinberg's approach [|l7|-pcfl 
suggests that renormalized coupling constants should be well-behaved. 

Below we use Weinberg's power counting with BPHZ renormalization procedure and compare the 
results of subtractively renormalized and cutoff EFT for 1 Sq NN scattering amplitude. We find that the 
results of two schemes coincide up to the order of the accuracy of given calculations. We also explore the 
possibility of non-perturbative renormalization via absorbing part of divergences into the two available 
"bare" coupling constants with successive removal of regularization. 

In this letter we quote the results, details of calculations will be presented elsewhere. 

We iterate subleading order potential substituting it into Lippman-Schwinger equation. To perform 
subleading order analysis we include two-derivative contact interaction part of the potential for 1 5o NN- 
scattering as it is enhanced in comparison with the two-pion exchange part || (note that to renormalize 
the amplitude obtained by iterating leading order potential no contact terms with derivatives are required 
H, JO]) and therefore the sub-leading order potential has the following form: 



V(p',p) = C + C 2 (p' 2 + p 2 ) + V4 (p',p) 



where C and C 2 are ("bare") coupling constants, and 



(p' -P) +ml' 
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We need to extract the 5-wave part from the solution to the Lippman-Schwinger equation for the TV TV 
scattering amplitude T: 



T(p',p) = V(p',p) 



where 



M 



d 3 q V (p',q) G(q) T(q,p) = V(p',p) + V (p',q) ® G(q)®T(q,p) , 



(3) 



G(q) = 



k A -q 2 



), M is the nucleon mass and p is on-mass-shell 



integration over "q - variable" J d 3 q is denoted by 

value of three-momenta. For physical amplitude p' = p. 

When substituted into Lippman-Schwinger equation (||), potential V (p',p) generates divergences. To 
regularize potential we introduce factorized form-factor into the contact interaction part of V (p' ,p) and 
write: 



V R {p',p) 



C + C 2 [p 



P 



x(A,p')x(A,p) + Vir(p',p), 



where x is a regularization function. For definiteness let us use 

A 4 



x(A,p) 



(A 2 + p 2 



(4) 
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Using Vr(j>',p) instead of V(p',p), we rewrite the Lippman-Schwinger equation (||) into the matrix 
form: 



T(p o; p) 'SH^>(K 

^/a{x(a.„<)(;-; 



c, c 2 
c 2 , 

c, c 2 
c 2 , 



1, 
p 2 , 

1, 
q 2 , 



x(A,p) 



X(A,<? 



vAp',p), 
0, 

V v (p',q), 







0, 
G(q), 0\ ( T(q,p), 



0, 



0, 



(0) 
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We solve (||) for the 1 Sq amplitude and obtain: 

T ( ) = T + [ g + 2 P 2g 2 + C 2 G 3 {p)] V^p) 2 - [1 + C 2 G 2 {p)] + GfG^)^) 

* [l + C 2 G 2 {p)f -[C + 2p 2 C 2 + C 2 G 3 {p))G l {p) 

where T n is 1 So wave solution of Lippman- Schwinger equation (||) for V (p',p) — 14 (p',p). Functions 
Gi(p) are given by infinite series as follows: 

Gi(p) = X 2 (A, q) ® G (g) + X (A, g x ) ® G ® K (g x , g 2 ) ® G (g 2 ) ® x (A, g 2 ) + 

X (A, 9l ) (8) G ® 14 (g 1( g 2 ) ® G (g 2 ) ® 14 (92,93) ® G (g 3 ) ® X (A, g 3 ) H (8) 

G 2 ( P ) = x 2 (A, g) ® 1 + x (A, gj ® 14 9 2 ) ® G (9a) ® X (A, 9a) + 
X (A, q x ) ® 14 (91, 9a) ® G (g 2 ) 14 (g 2 , g 3 ) ® G (g 3 ) ® x (A, g 3 ) + • • • (9) 

G 3 (p) = X 2 (A, g) ® G- 1 (g) + x (A, gj ® 14 (g l5 g 2 ) ® x (A, g 2 ) + 

X (A, 9l ) ® 14 (g X) g 2 ) ® G (g 2 ) ® V* (g 2 , g 3 ) ® x (A, g 3 ) + • • ■ (10) 

and the functions Vi(p), V 2 (p) satisfy integral equations: 

Vi(p) = x (A,p) + ^ (g) G (g) ® 14 (g, p) , (11) 

V 2 (p) = x (A, 9) ® 14 (9, P) + ^2 (g) ® G (g) ® 14 (9, P) • (12) 

Expression (|7|) for the amplitude remains valid in the framework of dimensional regularization provided 
that integrals are regulated using this regularization and X - function is substituted by 1. 

To give physical meaning to Eq. (7) we need to get rid off divergences, which show up in A — * 00 limit. 
Note that the amplitude can be made finite using non-perturbative renormalization, i.e. renormalizing 
two available "bare" parameters G and G 2 and taking A — > 00 limit in non-perturbative expression of 

Let us describe this procedure . We fix G and G 2 from two conditions of finitcness of the amplitude 
at p 2 — and p 2 = — /i 2 (p = ifi). As T^{p) does not contain divergences, we require the finiteness of 
a(p) = l/[T (p)-T n (p)}: 

rm [1 + G 2 G 2 (0)] 2 - [G + G 2 2 G 3 (0)] Gi(0) 

a[U) [C + G 2 G 3 (0)] V^O) 2 - 2C 2 V 1 (0)V 2 (0) [1 + G 2 G 2 (0)] + Gfd (0)1^(0) ' ( 6) 

, , = [1 + G 2 G 2 (z M )] 2 — [C — 2 M 2 G 2 + CjGsM] G^) 

a(N [C -2^ 2 C 2 + C 2 G 3 (i^]V 1 (i^ 2 -2C 2 V 1 (i^V 2 {tfi)[l + C 2 G 2 (in)} + C 2 G 1 (i^V 2 {ifi) [ ' 

In - @ o(0) and a(/i 2 ) are finite. 

Analyzing expressions ([§l)-(|l2|) it is straightforward to obtain following relations: 

Gx(p) = a A + 01 In A + a 2 (p) + (1/A) (15) 

G 2 {p) = b A 3 + hA 2 + b 2 (p)A + b 3 {p) In A + b 4 {p) + (1/A) (16) 

G 3 ( P ) = c A 5 + Cl A 4 + c 2 (p)A 3 + c 3 (p)A 2 + c 4 (p) In A + c 5 (p) + (1/A) (17) 

^i(p) = Vi(p) + 0(l/A) (18) 

F 2 (p) = wA + d(p) + (1/A) , (19) 

where ao,ai,bo,bi,ci,ci,v are constants independent of A and a 2 (p),b 2 (p),b 3 (p),c 2 (p) — 
cs{p), Vi(p), d(p) are some functions of p also independent of A. 



Substituting expansion (15)-([L9|) into expressions for a(0) and a(/i) we solve the system of two equations 
for G and G2 . Substituting these solutions for G and G2 into expression for the amplitude (Q) and taking 
into account expansions (|lq)-()l9|) we obtain in the limit A — > 00 amplitude which is finite: 



3 



T (p) =Ti I (p) + — (20) 
a(p) 

where 

a(p) = {(ai^Vfiti) + a 2 (p)) (c 2 (p) - c 2 (0) - 4b oP 2 ) - (a(O)lf(O) + a 2 (0)) x 

(c 2 (p) - c 2 ( M ) - 4& P 2 - 26 ^ 2 ) + 02 (p) (c 2 (0) - c 2 (/i) - 26 M 2 )}/ [(c 2 (/x) - c 2 (0) + 26 oA i 2 ) V?(p)] (21) 



When using the same renormalization condition and dimensional regularization we obtain finite re- 
sult which differs from (]2l]) in analogy with the case of pionless EFT [E2j. This difference between two 
results is due to the fact that dimensional regularization discards all power-law divergences, i.e. dimen- 
sional regularization is in fact not only regularization in standard sense but it also subtracts power-law 
divergences. 

As it was argued in pjJ, the above described renormalization scheme is not consistent with EFT. To 
carry out the consistent renormalization procedure it is necessary to absorb all divergences into coupling 
constants or in other words, subtract all divergent integrals. To do so let us start from the observation 
that G 2 , G 3 and V2 can be expressed as follows: 

G 2 ( P ) = A 2 + A 7I G 1 (p), 
G 3 ( P ) = A 3 + A T G 2 (p), 

V 2 (p) = A n V 1 (p), (22) 

where 

A 2 = x 2 (A,q)®l, A 3 -x 2 (A, < 7)®G- 1 ( < z), A w = % (A, q) ® V„ (q, p) . (23) 

Relations ( p2| ) are correct up to terms which are vanishing after renormalization is performed and regu- 
larization is removed , i.e. after sub-divergences and overall divergences are substracted. Note that these 
relations become exact in dimensional regularization provided that in (^3|) x~f ac tors are substituted by 
1. 

Using relations (]23) we express amplitude (R) as 



Tn=TM , Vi(p) [C - A„ G 2 (2 + A 2 G 2 ) + 2 p 2 G 2 + A 3 G 2 ] 

0[P> wW 1 + 2A 2 G 2 + A 2 G 2 - Gi(p) [G - A, G 2 (2 + A 2 G 2 ) + 2 p* G 2 + A 3 G|] 1 ' 

To renormalize amplitude (|24|) following standard BPHZ renormalization scheme we subtract divergent 
parts from G\(p), A 2 , A 3 and A n and substitute G and G 2 by their renormalized values Cr and C 2 r 
and obtain: 

T R ( ) = T ( ) + V?{P) ~ A - ° 2R (2 + F2C2R) + 2 ^ ° 2R + {F31 + f32 P ^ C ^ 

(P> v(P) 1 + 2F 2 C 2R + FiC 2 2R - G R {p) [C R - A R C 2R (2 + F 2 C 2R ) + 2 p* C 2R + (F 31 + F 32 p*) C\ R \ 

(25) 

Here F 2 , F 31 and F 32 are (arbitrary) finite parts. Note that arbitrary finite parts are also present in A^ 
and G R (p), the particular choice for the finite parts corresponds to the particular choice of renormalization 
condition. 

Expression ( [25] ) agrees with the one for the renormalized amplitude, obtained in the framework of 
dimensional regularization Q . Amplitude obtained in j^] can be reproduced from Eq. ( p5|) by appropriate 
choice of finite parts of G\(p) R and A R , and by the conditions F 2 = F 3 i ~ F 32 = 0. Hence, as it 
was expected, dimensional and cutoff regularizations lead to the same results provided that the same 
renormalization condition is used. 



Let us now turn to the cut-off EFT. In this approach we absorb part of divergences into the "bare" 
coupling constants G and G 2 but instead of considering A — > 00 limit, we choose the value of cutoff param- 
eter A in such a way that the amplitude of cut-off theory reproduces the above calculated subtractively 
renormalized amplitude up to sub-sub-leading order. 
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We determine C and C 2 from two conditions, the first of which is to absorbs all terms in Eq. (24) for 
p = 0, which are divergent if the A — * 00 limit is taken. This condition can be written as: 

1 + 2A 2 C 2 + A 2 2 C 2 2 - G^p) [C - A, C 2 (2 + A 2 C 2 ) + 2 p 2 C 2 + A 3 C 2 ] _ 



V 2 {Q) [C 



1 + 2A 2 rC 2 r 



A 2 r 2 



A, C 2 (2 + A 2 C 2 ) + 2p 2 C 2 
-Gf{p) [Cr-A^r C 2fl (2 4 



V A 3 C 2 2 ] 

A 2 rC 2 r) + A 3fl C| K ] 



V 2 (0) [C R - A^r C 2R (2 + A 2R C 2R ) + A 3R C. 



(26) 



2Ri 



Another condition can be obtained by expanding Eq. ( p4| ) in powers of C 2 and requiring the finiteness of 
the first order term in A — > 00 limit. Taking into account Eq. (J2q) , this requirement can be expressed as 



C2 
C 2 



C 2 R 



(27) 



Substituting the solutions to Eqs. (|26|)-(|27|) into the (unrenormalized) solution of Lippman-Schwinger 
equation given by Eq. (0) and expanding in powers of C 2 r we obtain: 



nip) = t„( P ) + - 

2V 1 (p) (A^ip) 



V 2 ( P ) Cj- 



2 VAp) (p 2 Vi(p) - A« Vi(p) - C R Af V x ip)) C 2R 



- C R (Ax - 
V 2 (p))C 2R 



Gib) ) (1 
2V 1 (p)C 2R Cr{A 2 V 1 ( P )- 



-C R [A 1 -G 1 {p)}) 
G 2 {p)V 1 {p)+G l {p)V 2 {p)} 



(I + CrIAx-G^pW 



Vi(o) (i + AxCr) (l + CflfA-dW) c R y 



(28) 



where Ai contains divergent part of G\ as well as finite part, the latter depending on the choice of 
renormalization scheme. Note that fourth and fifth terms of Eq. ( ^8|) are suppressed by an additional 
factor of 1 /A in comparison with third term. Higher order terms in C 2 r which are indicated by dots are 
suppressed by higher orders of small parameter (three-momentum p, pion mass or the renormalization 
scale which we choose of the order of p) . As our potential includes only leading and subleading order terms, 
calculated amplitude can be trusted up to this order, provided that formally higher order contributions 
are indeed suppressed. Note that these higher order terms also contain factors like p 4 A 2 C 2R . To suppress 
this kind of contributions as well as terms with inverse powers of A, we need to choose cutoff A of the 
order of large scale parameter Q which governs the behavior of Cr, C 2 r ||, |D|. With this choice of 
cutoff parameter one can easily see that two amplitudes given by Eq. (|2^) and Eq. ( |28| ) coincide up to 
the order of accuracy of our calculations the difference being of higher order. 



In this work we have iterated cutoff regulated 1 Sq partial wave NN potential of subleading order 
EFT and wrote a closed expression for the scattering amplitude. We have performed non-perturbative 
renormalization via absorbing part of divergences into two available coupling constants and taking the 
removed regularization limit. Within this renormalization scheme cutoff and dimensional regularizations 
both lead to finite results but these results differ from each other. This difference between two results is 
due to the fact that dimensional regularization discards all power-law divergences. We argue that this 
kind of renormalization is not consistent with EFT approach. 

We advocate the point of view that Weinberg's power counting should be applied to renormalizcd 
diagrams ||, 0, |L5| |, p6[, i.e. one should draw all (an infinite number of them for NN-scattering) 
relevant diagrams, renormalize and sum them up. This summation is performed by solving Lippman- 
Schwinger equation (or Schrodinger equation). In fact we have solved the regularized equation and 
performed renormalization by subtracting divergences in the solution. For the case considered in this work 
it is not difficult to demonstrate that expanding the renormalized solution to the Lippmann-Schwinger 
equation we reproduce separate renormalized Feynman diagrams. Of course in general the solutions to 
Lippmann-Schwinger equation (or Schrodinger equation) can contain more contributions than just the 
sum of (an infinite number of) perturbative diagrams, but note that there is no power counting for such 
non-perturbative contributions. 

We have implemented the renormalization procedure using BPHZ procedure. Within this scheme 
no counter-terms are included in the Lagrangian and therefore it is free of formal inconsistences that 
to remove divergences from diagrams of given order calculations one needs to include contributions of 
counter-terms which are of higher order. Let us emphasize that we have not attempted to solve this 
problem of formal inconsistences. We argue that this issue is irrelevant and it can be avoided by applying 
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power counting to renormalized diagrams. The usage of BPHZ renormalization proves to be very useful 
in dealing with this problem. 

The main result of our work reads as follows. While the power counting is consistently formulated in 
subtractively renormalised EFT, in most cases we are not able to solve regularized equations and therefore 
it is technically very difficult (or rather impossible) to implement subtractive renormalization scheme. 
In practice most of practically important calculations are performed using cutoff EFT. For the potential 
considered in this work we have compared the amplitude of cutoff EFT with subtractively renormalized 
amplitude. The cutoff amplitude is equal to subtractively renormalized amplitude up to (including) sub- 
leading order, i.e. the difference between these two is of higher order than the order of the accuracy of 
our calculations (note that subtractively renormalized subleading order calculations reproduce the phase 
shift analyses for 1 Sq NN scattering reasonably well ]l4|]). 

We consider the present work as a demonstration of a general feature of cutoff EFT that it reproduces 
the results of subtractively renormalized theory up to the order of the accuracy of given calculations. 
This suggests that the considerable success of Weinberg's approach implemented in cutoff EFT 0, ]Tq| , 
Jlj| , |20fj is by far not accidental. 
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